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Abstract 

We consider the polar factorization of vector valued mappings, introduced in [Sj, 
in the case of a family of mappings depending on a parameter. We investigate the 
regularity with respect to this parameter of the terms of the polar factorization by 
constructing some a priori bounds. To do so, we consider the linearization of the 
associated Monge- Ampere equation. 

1 Introduction 

Polar factorization and Monge-Ampere equation 

Brenier in jHj showed that given Q a bounded open set of R d such that \dQ\ = 0, with |.| the 
Lebesgue measure of IR d , every Lebesgue measurable mapping X e L 2 (f2,IR d ) satisfying 
the non- degeneracy condition 

(1) Vfi C R d measurable, \B\ = =^ |X -1 (£)| = 
can be factorized in the following (unique) way: 

(2) X = V$og, 

where $ is a convex function and g belongs to G(Q) the set of Lebesgue-measure pre- 
serving mappings of Q, defined by 



(3) g e G(Q) V/ G C b (Q), / f(g(x)) dx = / f(x) dx, 

Jn Jn 

where Cf, is the set of bounded continuous functions. If da denotes the Lebesgue measure 
of Q, the push-forward of da by X, that we denote X#da, is the measure p defined by 

(4) VfeC b (R d ), [ fdp= [ f(X(a))da. 
jR d Jn 
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One sees first that the condition (JTJ) is equivalent to the fact that p is absolutely continuous 
with respect to the Lebesgue measure, or has a density in L 1 (M c( , dx). Then <fr satisfies in 
Q the Monge- Ampere equation: 

p(V&(x))detD 2 <f>(x) = 1 

in the following weak sense: 

(5) WgeC b (R d ), [ g{V*(y))dy = f g(x)dp{x). 

Jn Jm d 

VP, the Legendre transform of defined by 

(6) *(y) = sup{x • y - $(ar)}, 

x&n 

satisfies the Monge- Ampere equation 

det D 2 &(x) = p(x) 

in the following weak sense: 

(7) V/ G C 6 (n), / /(V*(x))dp(x) = / f(y)dy. 

JK. d Jn 

Note that the existence and uniqueness of the pair V<&, V\l/ and the validity of (JHJ) is 
not subject to the condition (JTJ) (see |22| Th 2.12 for this precise fact, and for a complete 
reference on polar factorization and optimal transportation). However may not hold. 
Note also that this formulation of the second boundary value problem for the Monge- 
Ampere equation is strictly weaker than the Aleksandrov formulation (see jH] where the 
different formulations are compared and where it is shown that they may not coincide if 
some extra conditions are not satisfied). 

The periodic case The polar factorization of maps on general Riemannian manifolds 
has been treated by J7j, and also in the particular case of the flat torus by JU]- Given 
X a mapping of T d = M. d /Z d into itself, we look for a pair ($, g) such that 

1. g is measure preserving from T d into itself, 

2. $ is convex from IR d to K and $ — \x\ 2 /2 is periodic, 

3. X = V3> o g (Note that the condition above ensures that V<fr — x is Z d periodic). 
Then under the non-degeneracy condition {IJ), there exists a unique such pair (g, V«&). 
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Introducing the time-dependence 

In this paper we are interested in the following problem: given a "time" dependent family 
of mappings t — > X(t, .), where for all t, X(£) maps fl in IR d , we investigate the regularity 
of the curve t -> (g(t, .), *(t, .), *(*, •))• 

We state different results under different assumptions. The weakest assumption is that 
p = 'K^da, X and <9 t X belong to L°° in time and space. In this case 9 t V$ and <9 t g are 
bounded as measures (Th. 12. 

Under the additional assumption that p is close to 1 (or actually to a continuous 
positive function) in L°° norm (but we do not ask for continuity), we obtain that dt<& 
belongs to C a for some a > (Th. 12. 2j) . To this purpose we use a local maximum principle 
for solutions of degenerate elliptic equations (Theorem 13.51 Theorem I3.7J1 obtained by 
Murthy and Stampacchia (|18j) and Trudinger (|20j). and use a result by Caffarelli and 
Gutierrez (jHj) that establishes the Harnack inequality for solutions of the homogeneous 
linearized Monge- Ampere equation (Theorem 13.4)1 . 

The polar factorization has the following geometrical interpretation: if X = V<I?og, as 
in (J2J), then g is the projection, in the L 2 (fl, M d ) sense, of X on G(fl), the set of Lebesgue 
measure preserving mappings. Therefore our study amounts to examine the continuity 
and the differentiability of the projection operator on G(Q). We also briefly discuss a 
variant of the Hodge decomposition of vector fields that appears naturally in this study. 

Our results have an immediate application to the semi-geostrophic equations, a system 
arising in meteorology to model frontogenesis (see J2])- They allow in particular to define 
the velocity in the physical space, a fact that was not known for weak solutions. We discuss 
this application in a more extensive way in sectional 



1.1 Heuristics 

We present here some formal computations, assuming that all the terms considered are 
smooth enough. Suppose that fl is bounded, and for any t we denote by dp(t, •) = 
X(t, -)#G?a (with da the Lebesgue measure on fl) the measure defined by (J3J). Then for 
all t, *(t, •), *(*, ■) are as in Ipf7jl . 

Parallel with the Hodge decomposition of vector fields 

By differentiating (j2J) with respect to time one finds 

d t X(t, a) = dtV*(t, g(t, a)) + D 2 <f>(t, g(t, a))d t g(t, a). 
If X is invertible, one can write 

(8) d t X.(t,a) =v(t,X(t,a)) 

for some "Eulerian" vector field v(t, x) defined dp a.e. Note that p = X^rfa and v will 
be linked through the mass conservation constraint 

(9) d t p + V ■ (pv) = 0. 
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g will then also be invertible and composing with g 1 one gets: 

(10) v(t, V$(i, x)) = $V*(t, x) + D 2 ${t, x)w{t, x) 

with w = d t g(t, g _1 (t, x)). Since for all t, g(i) G G(fi), it follows that w is divergence 
free. Composing with = V3> 1 we obtain 

v = <%V$(V*) + D 2 ® w(V*). 

It is easily checked that w = D 2 <f> u>(V\E') satisfies 

V • (pw) = 0, 

therefore the second term in the decomposition (fTUj) does not move mass. It plays the 
role of a divergence free vector field for a uniform density. 

Note that a similar decomposition is performed in the study of the incompressible inho- 
mogeneous Navier-Stokes equation in ^1 where for a given velocity field v , and a density 
p > 0, one seeks to decompose v as 

v = — Vp + w, V • w = 0. 
P 

The next proposition shows that, in the non-degenerate case where $ is smooth and 
strictly convex, the decomposition (jlOj) is defined in an unique way. 

Proposition 1.1 Let v G L 2 (R d } dp; M. d ), let $ : — > R d be C 2 and strictly convex on f2 ; 
with p = V$#c/a. Then there exists a unique decomposition of v such that 

(11) u(V*) = Vp + D 2 <&w 
with (Vp, w) G L 2 (fi; R d ) ; V • w = 0, w ■ dtt = 0. 



Proof: We only sketch the proof of this classical result, w can be found by looking for 



(12) 



inf 

w e L 2 {Q;1 

w-dn = o 

V • w = 



</5 



to* ■ D 2 <& ■ w - u(V<&) • to}. 



Using the strict convexity of $ we have D 2 <I> > A/ on Q, and we obtain that 



\w\ 



L 2 {Q) 



2 

< - 
- A 



p\v\ 



1/2 



The functional to minimize is strictly convex, and weakly lower semi continuous, therefore 
the problem admits a minimizer. For the uniqueness of the decomposition, notice that if 

= Vp + D 2 ®w 

for Vp, w G L 2 , multiplying by w and integrating over Q, we get that Vp, w = 0. There- 
fore, if v governs the evolution of p through the equation Q, the decomposition (fTTj) will 
coincide with flTUj) and will yield Vp = d t V$>. 
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The associated elliptic problems: The linearized Monge- Ampere equation 

Multiplying (fTU|) by D 2 ^^ 1 , we find that d t & will be solution of the following elliptic 
problem: 

V • (D 2 * _1 VSi#) = V • (£> 2 *-y V*)). 

On the other hand, \I> = $* (see (|Sjl ) solves formally the equation 

det D 2 V = p. 

Then for any (d x d) matrices A, B we have 

det (A + tB) = det A + t trace(*A*5) + o(t) 

where A* is the matrix of cofactors (or co-matrix) of A and thus, formally, d t ^ solves the 
elliptic equation 

Mijdijdt® = d t p, 
where (My)^^..^ is the co-matrix of D 2 *&, given by 

M = detD 2 *^ 2 *]" 1 = p£> 2 $(V*). 

Then if M is the co-matrix of a second derivative matrix, for all j G [l..d] 

d 

Y^diMaix) = 0, 

i=l 

and using this and the equation Q, we obtain a divergence formulation of the problem: 
(13) V -(MVd t V) = d t p = -V -(pv). 

In the case where p is smooth and supported in a convex set, it will be shown using classical 
elliptic regularity and results on Monge- Ampere equation, that the decomposition holds 
(Proposition I4.1J1 and that the terms are smooth. 

For a generic, non-necessarily smooth p, we see that the difficulty will be coming from the 
lack of regularity and ellipticity of this equation. Indeed we only know a-priori that D 2 $ 
is a measure. If p is close to 1 in L°° norm, we get that D 2 & is in l7 loc for some p < oo, 
and thus non necessarily uniformly elliptic. 

2 Results 

Notations 

In the remainder of the paper Q will be kept fixed once for all and chosen bounded and 
convex. We will furthermore assume for simplicity (although one may possibly remove 
this assumption through approximation) that it is smooth and strictly convex. 
The Lebesgue measure of Q, Xn£ d , will be denoted in short da. 

For compatibility p will be a probability measure on M. d and Q of Lebesgue measure one. 
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M.(Q) will design the set of (possibly vector valued) bounded measures on Q, with norm 

ll-ll-M(ft)- 

For M a (d x d) matrix, and u, v two vectors of M. d , uMv will denote £\ ■ UiM^Vy 
I will be an non-empty open interval of R. 

We still use dp(t, •) = X(£, -)#da, the functions 3>(t, •), &(t, •) will be as in (0 Ej) with 
(p(i, Since they are defined only up to a constant, we will impose the condition: 

(14) Vt 6 I, / *(t,x) dx = 0, 

Jn 

and this sets also ^ through the relation ^ = 



Theorem 2.1 Lei fi, / be as above, let X : / x fi — > R d . Let, /or any t & I, dp(t, •) = 
X(t, ')jfda as in ^jj. Assume that (X, c^X) G L°°(/ x fi), wift -R = ||X||ioo( JX Q) ; and 
assume that p G L°°(I x R d ). Ta/ce 

X(*) = V*(t) o g(t), g(t) = V*(t) o X(t) 

to 6e i/ie polar factorization o/X as m (0) where we impose \14\) . Then 

1. for a.e. t G /, ^V$(t, •) is a bounded measure in Q with 

i 

||9tV$||i«>(/ i> M(r2)) < C(-R, rf, fi)||p||| / oo( JXjBiJ )||<9 t X|| i oo (/XjBr ) 

andd t & G L°° (I , L 1 * (tt)) with 1* =d/(d-l). 

2. $ fresp. belongs to C a (T,C°(n)) (resp. to C a (I; C°(B R ))) for some a G]0,1[. 
5. For a.e. £ G I, <9tg zs a bounded measure on Q with 

||^g|U°°(/;M(n)) < C(-R, rf, fi)||p||ioo( /XjBfl )||(9 t X||z / oo( /X Q). 

^. L/p is supported in Cl' for some open set Q' , and < A < p(-, ■) < A on Q' , for some 
(A, A) G R+, then there exists j3 g]0, 1[ such that for any uj' CC Q! , 

V* G C^I; C°(u/)), 

with 13 depending on A/A. 
5. If in addition Q' is convex, then there exists f3' g]0, 1[ such that for any w CC H, 

V$ G C p '{I-C°{oo)). 



Theorem 2.2 Under the assumptions of Theorem \2.1\ and assuming that p is supported 
in Q' , for some open set Q' , we have: 
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1. There exists e > such that if \p — 1| < e < e in Q! , then there exists a > 
(depending on e) such that, for any w' CC Q' , 

d t Ve L°°(I; C a (uj')). 

If in addition Q' is convex, for any w CC Q, 

d t $ G L°°(I; C a {ijj)). 

2. For any p < 2, there exists e(p) > such that, if \p — 1| < e(p) in Q' , for any 

w' cc n', 

d t V* G L°°(J; L p (uj')). 
If in addition Q' is convex, for any w CC Q, 

d t V$ E L°°(/;L p (cj)). 



Remark: The Theorem remains true if one replaces the condition \p— 1| < eby \p — f\ < e 
with / a positive continuous function and the bounds will then depend on the modulus 
of continuity of / (see jl]) . 



We also state the result in the periodic case: In this setting we have the following 
theorem, which is just an adaptation of the two previous: 

Theorem 2.3 Under the assumptions that p 6 L°°(I x T d ), d t X G L°°(I x T d ), we have: 

1. With the same bounds as in Theorem \2.1\ 

d t V$ G L°°(/;7W(T d )), 
dtS G L°°(T, A4(T d )), 

and for some a > 0, we have 

G C a {I; C°{T d )). 

2. If for all (t, x) E (I x T d ) we have < A < p(t, x) < A ; then for some (3 > 
depending on (A, A) G M+, 

g,V$,V* G C^(T,L°°(T d )). 

3. There exists eo such that if \p — 1| < e < eo, then for some a > depending on e, 



dtV G L°°(J; C a (T )) 



9t* G C a (T d )). 
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4- For any p < 2 there exists e(p) such that if \p — 1| < e(p) then 

d t W E L°°(J; L p (T d )), 
d t V® e L°°{I- L p (T d )), 
d tg eL™(I;L p (T d )). 

Remark: in this case, the absence of boundary allows to have a bound over T d and not 
only interior estimates as in the previous results. 



2.1 Related results 

The linearized Monge-Ampere equation 

The linearized Monge-Ampere equation (LMA) is a well known equation, since it is used 
to carry out the continuity method, in order to obtain classical solutions of the Monge- 
Ampere equation (see [13], chapter 17). However for this purpose this is always made in 
the case where the densities and the domains considered are smooth, and thus the LMA 
equation is uniformly elliptic. 

In the non-smooth case, [§] proved Harnack inequality for solutions of 

MydijU = 

with M the co-matrix of D 2 ^, for some \l/ convex, under the assumption that the measure 
p = det D 2 ^> satisfies the following absolute continuity condition: 

C: For any < 5\ < 1 there exists < 62 < 1 such that for any section S and any 
measurable set E C S, 

(15) if ||f < S 2 then 4|r < *i, 

\S\ p(S) 

(a section is a set of the form 

S t (x ) = {x\^(x) - *(x ) <P'(x-x )+t, pe d^(x )}). 

They showed that the solution of (det D 2 ^/)(D 2 ^)~j 1 D i jU = satisfies a Harnack in- 
equality on the sections of ^ and subsequently is C a . The precise result is stated below 
(Theorem 1.14)1 . We will use this result to obtain the first part of Theorem 12.21 Note that 
the condition implies C 1,a regularity of the Aleksandrov solution of det D 2 ^> = p 
([S3)- Note also that the condition (fl3|) is satisfied when the density p is bounded between 
two positive constants. We will also obtain some results (Theorem 12. lj) in the degenerate 
case when the condition ([T5|) is not satisfied and show in some counterexamples (section 
IE} that when this condition is not fulfilled, the result of Theorem 12.21 does not hold. 

Maximum principles for degenerate elliptic equations 

We will use a local maximum principle for degenerate elliptic equations to obtain Holder 
continuity in Theorem 12.21 Consider the problem 



V • (M(x)Vu(x)) = V • f(x) 
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where M(x) = M^-(x), G [l..d] is a symmetric positive semi-definite, matrix, f(x) = 
(fi(x))i G [l..d\. In the cases we will study, we will not have the usual uniform ellipticity 
condition 



with / the d x d identity matrix, and for some positive numbers A, A, but a condition of 
the form 



for some non negative measurable functions X(x),A(x). Under the assumption that 
(A _1 (x), A(x)) G L p loc [yt) for some p > d and that / G L°°, we can obtain a bound 
on the solution u in L^ c . Properly localized, this bound with the Harnack inequality 
( Theorem 13. 4|) will yield Holder continuity of the solution of the LMA equation (JT3"j) . 
This type of maximum principles have been already obtained in [TH], [201, (see also [H]), 
and we will use them under the forms of Theorems 13. 5| 13 .7\ and Corollary 13.61 Note 
however that the condition (fi7)J) is not know by itself to guaranty Holder continuity of the 
solution, but only a L°° bound. 

It can be interesting to point out that we will thus use both the divergence and non- 
divergence structure of the LMA to obtain our results. 

3 Some preliminary results 

In this section we state the results that we are going to need for the proofs of the theorems. 
The reader may skip this section and come back to it whenever needed. Note that all 
these results can be extended to the periodic case. 

3.1 Regularity for solutions of Monge- Ampere equation 

Theorem 3.1 Let be bounded, C°° , strictly convex, and |0| = 1. Let p be a prob- 

ability measure in fl' , belong to C 00 ^'), and satisfy < A < p(t,x) < A for some pair 
(A, A). Then there exists a unique (up to a constant) solution of 

detD 2 * = p, 
V* maps Q 1 to Q, 

in the sense of F^j. The solution ^ belongs to C°°{Q!), and defined as in belongs 
toC°°{n). 

For this the reader can refer to jlJ-JH], [TB] . |2T| . 
The next Theorem can be found in j^J, [S], [7j. 

Theorem 3.2 Let p be supported in fl f with Q' open, satisfy < A < p < A, and let 
be solution of 



XI < M < AI 



(16) 



X(x)I <M< A(x)I 



det£> 2 * = p, 
V* maps Q' to Q 
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in the sense of TO with Q convex. Then for some a G]0,1[ depending on A/A, \E r G 
C^ip!). If moreover Vl' is also convex then (resp. its Legendre transform is in 
C 1 ' 01 ^') (resp. in C l > a (£l)). 



The next Theorem can be found in 

Theorem 3.3 Let Q be normalized so that B\ C Q C B^. Let \I/ be a convex Aleksandrov 
solution of 

det D 2 ^ = p, 
* = on on. 

Then for every p < oo there exists e(p) such that if \p — 1| < e(p) then G W^f^) and 

11*11 W*(B 1/3 ) < C(e). 

Remark 1: This implies also, maybe for a smaller value of e(p) that one can also have 
P 2 *- 1 ^^) < C"(e). 

Remark 2: The theorem remains true if one replaces \p — 1| < e by |p — f\ < e, for some 
continuous positive /, and the bounds depends on the modulus of continuity of /. 



3.2 The linearized Monge- Ampere equation 

We state here the result of jU] evoked in the previous section: 

Theorem 3.4 Let Q be a domain in ~R d , let U be an Aleksandrov solution in Q of 

det D 2 U = [i 

where [i the satisfies the condition hlfy) . Let w be a solution in Q of the linearized homo- 
geneous Monge- Ampere equation 

Aijdijw = 

where Aij is the co-matrix of D 2 U, let R > and y G Q be such that Br(jj) C f2, then 
for some (3 < 1 depending only on the condition < f73j) . for any r < R/4, 

osc(r/2) < (3osc(r), 

where 

osc(r) = M(r) — m(r), 

M(r) = sup w, m(r) = inf w. 

B r (y) Mv) 
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3.3 Maximum principle for degenerate elliptic equations 

We give here some results concerning degenerate elliptic equations of the form 
(17) V • (M(x)Vu(x)) = V • f(x) 

where M is symmetric non-negative matrix, / = i = l..d. The equation can be writ- 
ten di(Mijdju) = difi with summation over repeated indices. The usual strict ellipticity 
condition 

A|e| 2 < A%&£; < Aiei 2 for all f G R d , 
is replaced by the following 

d 

\ M ij\ + \ Mij \ e L lc( Q ) for some P, 

where M u denotes the inverse matrix of M. This is equivalent to the condition that there 
exists X(x),A(x) such that A -1 , A are in Lf oc (f2) and such that X(x)I < M(x) < A(x)L, 
in the sense of symmetric matrices. 

The class of admissible test functions is 

C(0) = {ve Wo'\n), M 1/2 Vv G L 2 {Q)}. 

A subsolution (resp. supersolution) u of (fTTj) is defined by the condition that for all 
non- negative v G C(Q), 

VvMVu - Vv ■ f < (>)0. 
Then, following [THj and we have the following results: 

Bound for Dirichlet boundary data 

We denote by S£ the set of d x d non negative symmetric matrices. 

Theorem 3.5 Let M : Q —> be such that M~ l is in L p {Vt\ Sj) for some p > d. Let f 
be in L°°(f2; Let u be a subsolution (supersolution) of 

V ■ (M(x)Vu(x)) = V- f(x) 

in Q, satisfying u < (u > 0) on dQ. Then 

supu(-u) < C(\\u + (u')\\ L a 0{n) + ||/|| £ oo(n)) 
n 

where C, C depends on a > 0,p > d, ||M _1 || £ p(n). 



This maximum principle can be precised in the following corollary, that will be crucial for 
the proof of Holder continuity in Theorem 12.21 

Corollary 3.6 Under the previous assumptions, for y G ; Br(v) C Q, if u is a subso- 
lution (supersolution) in Br of \17\) and u < (u >0) on 8Br, then 

supu(-u) < C\\M- l \\ LP{BR) \\f\\ L ^ [BR) R 5 , 

B R 

where 5 = 1 . 

P 
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Bound without boundary data 

Here we state a maximum principle that does not depend on the boundary data. Note 
that here we need to control the norm of both M and M~ x whereas we only needed to 
control M _1 above. 

Theorem 3.7 Let M : Q -> S% be such that M, M" 1 are both in L^ oc (Q), with p > d. 
Let f be in L°°(Q). Let u be a subsolution of 

V • (M(x)Vu(x)) = V • (/(:r)) 

in Q. Then we have for any ball B 2 r CC O and a > 

SUp U < Ci\\u + \\L*a(B 2R (y)) + C 2 k 
B R {y) 

where k = ||/||x,°°(b 2H ), Ci,C 2 depend on R,a ,p, ||M|| LP(B2H ), ||M _1 || LP ( B2fl ). 

3.4 Convex functions and Legendre transforms 

We state first the following classical lemma on convex functions: 

Lemma 3.8 Let if be a convex function from ~R d to R, globally Lipschitz with Lipschitz 
constant L. Then we have 



\D'M\m(b r ) < C(d)R d - 1 L 



Proof: we have 



\DMm{b r) < C [ A V 

J Br 

= / Vip ■ n 

JdB R 

< C(d)R d - 1 L. 



□ 



We recall here some useful properties of the Legendre transform. Let H be a convex 
domain, let : i— > R be C 1 convex. Let 0* be its Legendre transform defined by 

4>*(y) = supx -y - (j)(x). 

Then, for all x G f2, 

V0*(V0(a;)) =x. 

If moreover is C 2 strictly convex, then, for all x E VI, 
(18) D 2 (j)*(V<i)(x)) = D 2 (f)- 1 (x). 



From this we deduce the following lemma: 
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Lemma 3.9 Let Q be convex, let (t,x) i— > <fr(t,x) : / x Q i— ► R and (t,y) \—> *S?(t,y) : 

I xR d ^R be such that 

1. V$ (resp. V^f) belongs to C l (I x Q) (resp. belongs to C l (I x R d ) ), 

2. for all t E I, <&(t, •) is convex and Vl/(t, •) is the Legendre transform of <&(i, ■). 
then for every (t, x) 6 I X fi, 

(19) x) + *(t, V*(t, z)) = x • V$(t, x), 

(20) + 9i*(V$) = 0, 

(21) <9tV$ + £> 2 $<9 t V*(V$) = 0. 

Proof: the first identity expresses just the fact that $(t, ■), Sfr(t, •) are Legendre transforms 
of each other (see (jSJ)), then the two other come by differentiating with respect to time 
and then to space. 

□ 



4 Approximation by smooth functions 
4.1 Construction of smooth solutions. 

In this section we build an adequate smooth approximation of the problem. More pre- 
cisely, given a mapping X(£) and p(t) = X(t)#da, we construct an associated pair (p,v) 
satisfying 

(22) d tP + V • (pv) = 

and then find a "good" regularization of (p, v ). One of the problems is the following: it is 
known from a counterexample by Caffarelli (see |Hj), that when transporting a (smooth) 
density p\ onto another (smooth) density p2 by the gradient of a convex function, one can 
not expect the convex function to be C 1 unless p2 is supported and positive in a convex 
set. Therefore it is not enough to only regularize (by convolution for example) the density 
p = ~K#da, we must also approximate it by a density supported in a convex set. 
The density p and dtp are constructed from X, <9(X respectively by the following procedure: 

WfeCl(R d ), [ p(t,x)f(x)dx= [ f(X(t,a))da 
Jm. d Jn 

[ d t p(t,x)f{x)dx= [ V/(X(t,a)) ■ d t X(t,a)da. 
To define v such that d t p + V • (pv) = 0, we define the product pv as follows: 
V0 G Cf (/ x R d ; R d ), I pv -(pdtdx = [ 0(X(t, a)) ■ d t X(t, a) dtda. 

JlxR d JlxQ 

Since d t X e L°°, v is well defined dp a.e. and we have 

\\ v (ti -)IU <x '(R d ,<2 1 o(*)) — l|dtX(t, .)||i<x.(Q). 
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Now we construct (p n , v n ) a smooth approximating sequence for (p, v) as follows: (remem- 
ber that we have taken p(t, ■) to be supported in B R at any time t G I). We take 77 G C£° 
a standard convolution kernel, of integral 1, supported in B(0, 1) and positive. Take 
r] n = n d rj(nx). We also note xn+i/n the characteristic function of the ball B(0, R + 1/n). 
Let 

A 

Pn = {-XR+l/n +Vn* P)C n , 

n 

_ Vn * jpv) 
^n C n , 
Pn 

with c n chosen such that p n remains a probability measure. (Note that c n is close to 1 for 
n large). The purpose of this construction is to have the following properties: 

1- \\Pn,Vn\\L°° < \\p,v\\ L oo, 

2. p n ,v n satisfy the continuity equation (J22J), 

3. p n is supported and strictly positive in B(0, R + 1/n), and belongs to C°°(B(0, R + 
1/n)). 

4. If $ ra (t),^ r ra (t) are associated to p n (t) through (|5|7jl . then, for every t G /, $ n (t) 
converges uniformly on compact sets of Q to and *& n (t) converges uniformly 
on compact sets of M d to *&(t). This last result can be found in Therefore, 
9 t $ n , will converge in the distribution sense to d t &, d t ^f. 



Now we have the following regularity result, for smooth densities. Note that this result 
will only be used to legitimate the forthcoming computations, and not as an a-priori 
bound. 

Proposition 4.1 let I,Q be as above, let Q' be C°° strictly convex. For any t G /, let 
p(t, •) be a probability density in Q' , strictly positive in Q' with p G C°°(J x Q'). Let, for 
all t, <&(t, ■), •) be as in with (p(t),Q). Then, for any < a < 1, 

d t <S> G L°°{I, C 2 ' a {Q)), d t V G L°°(I, C 2 > a {tt')). 



Proof of Proposition Theorem 13.11 implies that for all t, D 2 ^ (resp. D 2 &) belongs 
to C°°(0') (resp. belongs to C°°(0)). 

Now we wish to solve det D 2i $?(t) = p(t) with t near to- We write a priori S&(t) = 
Sfr(to) + if — t )u + o(\t — to\), for some u, then we have 

det D 2 V(t) = det D 2 ^(t ) + (t - t )trace(MD 2 u) + o(\t - t \) 

where M is the comatrix of D 2 ^ defined by 

M(t, x) = det D 2 V(t, x) (D 2 V(t, x)) ^ . 
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Note that M belongs to C°°(Q') and is uniformly elliptic. Let us now show that d t ^ can 
indeed be sought as the solution of 

trace(M£> 2 w) = d t p 

with a suitable boundary condition. For this we introduce h a defining function for ft, 
(i.e. h G C°°(ft) is strictly convex and vanishes on <9ft, we can also impose |V/i|gn = 1). 
The condition VSt maps ft' on ft can be replaced by /i(V^) = on <9ft'. Now consider 
the operator 

T: 4> ^ (detD 2 iP,h(Vil;)\ d n>) 

defined on {ip G C 2,a (Q'),if) convex } and ranging in C a (ft') x C 1,a (dQ'). First note that 
a smooth solution of 



(23) *W = (p(*),0) 

will satisfy (JTj) and thus coincide (up to a constant) with *&(t). We now solve (J2~3*j) around 
to by the implicit function Theorem. The derivative of T at ^ is defined by 

dT(^)u = (I(u),B(u)) = (MijdijUihiiVtydiu) . 

The operator / = Mijdij is uniformly elliptic with coefficients My in C°°(f2 / ). We need 
also to show that the boundary operator B is strictly oblique: First, note that V/i = rii 
on dQ, where n\ is the outer unit normal to dfl. Moreover, if n,2 is the outer unit normal 
to dfl', it has been established in [S], [IS], EI], that there exists a constant C depending 
on Q, \\p\\c 2 (n'), an d therefore uniform on J, such that 

ft 2 -ni(V*) > C > 0. 

Thus the boundary condition is strictly oblique, uniformly with respect to t. It has been 
established in [T3j, p. 448, that the equation 

dF(V)u = (/i,0) 

with /i G C a (0') is solvable up to an additive constant if j Q , fi = 0. This condition is met 
by dtp, since f p(t,x) dx = 1. 

We conclude that the operator c?JF(^) is invertible on the set 



{{peC a (n'),Jp = 0}x{i 



0} 



ie. for each p G C a (0'), with f n ,p = 0, there exists a unique up to a constant solution 
it of dJ^iV^u = (/i, 0). Moreover, following [T3j, Theorem 6.30, it belongs to C 2,a {Q!). 
Therefore we can apply the implicit function Theorem and solve T(^(t)) = (p(t), 0) for t 
near t . By uniqueness of the solution of (JTJ), this solution will coincide with the solution 
of Theorem 13. II As we have built it, dt^(t, ■) — u is the unique (up to a constant) solution 
of 

(24) trace (MD 2 u) = d t p in Q', 

(25) Vw-nl(V*) =0 in 9ft', 
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and since d t p G C°°(Q'), d t *& belongs to C 2 ' a (Q') for any a < 1. 
We also have, using the identity (j2~Uj) 

d t $ + <9 4 *(V*) = 0. 

therefore d t $ G C 2 ' a (Cl) for any a < 1. 
This achieves the proof of Proposition 14.11 



□ 



5 Proof of Theorem 12.11 

Theorem 12. II will be deduced through approximation from the following proposition: 



Proposition 5.1 Let p satisfy the assumptions of Proposition [^77| above, with Cl' = Br, 
and ^ be as in 0). Let v(t,x) G be a smooth vector field on Br and satisfy on 
IxB R 

(26) d t p + V • (pv) = 0. 

Take 1 < p, r < oo, i + ^ — 1, q — Then for any t G I , for any w C O we have: 

(27) ||<9*v#||m ( uo < (IIpI^HI^'P^IU-^P^IUph) 172 , 



which implies in particular 

(28) ||atV*|Ui ( n) < C(R,d,Q) (||pM 2 ||l~) 1/2 

and for any t G I, for any to' C B R we have: 



(29) 



p\d t W\ q 



V? 



< ||p| V | 2 ||^p 2 *||^(^ 



1/2 



which implies in particular 



(30) 



p\d t v*\<c(R,dM\p\\I~ ( nM v \\\ 



2 II 2 



Proof of Proposition I5.il 

Using Proposition 14.11 we can perform the following computations. We have from (jSj) 

<9 t *p = / d t *(V$) 

Then we use the continuity equation: 

d t p + V • (pv) = 
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which implies for any smooth / 

/ fd tP = I pv V/. 

JR d 

We obtain 

' dt^dtp = I <9 4 V* • pv 

R d 

f. 

n 



dtV*(V*) ■ d t V<$> 
- / 9tV**(V$) ■ £> 2 $ ■ a t v*(v*) 



n 



where we have used (121 J) . Since we can write because this is a positive symmet 

matrix, we have 



JR d 



<9 4 V*(V$) .u(V#) 

= - / V£> 2 $a t v*(v*) • Vd 2 * _1 v(v*). 

This implies that 

(31) || v / d 2 ¥ a t v*(v$)|U 2(n) < iivWFVv*)!!^)- 

In order to estimate the right hand side, we write 

[ pv 1 ■ (D 2 <S>(W))~ l ■ v j 

jR d ) 

(32) < (||D 2 *|| L ,. (Bfl) ||pt; 2 ||^ (Bj?) ' 1/2 
In the second line we have used D 2 $(V*) = (D 2 ^)' 1 . From (J3TJ), 



1/2 



1/2 



llv 7 ^ ^v*!!^^ = ||v / ^$^v*(v$)|| i2(Q) 

< ||v^ 5 ¥" 1 v(v*)|| £a(n) . 



Writing 
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and, using Holder's inequality, we obtain for to C 

||atV*|U«( W) < ||V^¥ _1 a t v*|U 2(a;) ||^¥|U S (a;) 

< (\\p\v\%r iBR) \\D 2 n L r' { BjD 2 &\ 



1/2 



with g = 2T - . By taking p := s/2 we have 



||^V$|| L9H < ^\\p\v\ 2 \\ L r {BR) \\D 2 ^\\ Lr/{BR) \\D 2 ^\\ LP{uj) 
and g = This proves (|27)1. To obtain a bound on we write 



1/2 



'.D 2 $(V*)6> t V# 



/" a t v'*(v$) • £> 2 $ ■ ^v*(v$) 

in 



< WD^Wrr, ll " 1 - 121 



U r (B«)llPPl IIl'-'(Bh) 

from (jSI} and (jH2]) - Then using Holder's inequality, with q = ^f-, we obtain for a/ C -Br 



p|^V*| 9 



1/9 



< 



D 2 *(V*)8tV* 



1/2 



p|[L> 2 $(V*)] 



_l| S /2 



1/. 



The first factor of the right hand product has been estimated above, and the second is 
equal to I / p\D 2 ^f\ s/2 ) . We conclude that 



p\d t VV\ q 



< 



\D 2 V\W { B R )\\p\v\ 2 \\ L rl {BR) 



,1/2 



P 



D 2 V\ s/2 



l/s 



Taking again p := s/2, we have proved 
The bounds (1281 IHUj) are obtained as follows: we know from Lemma 13.81 that 

\\d 2 v\\ LHBr) <c(R,d,n), 

\\D 2 $\\ L i m <C{R,d,n). 

Taking in (|27l l2Uj) r = +oo,r' = l,p — 1 we obtain the desired bounds. This ends the 
proof of Proposition 15.11 

□ 



5.1 Proof of Theorem 12.11 

Proof of the bound on d t V& 

Here we prove points 1,2,4,5 of Theorem 12.11 To obtain point 1, we just need to pass to 
the limit in the estimate We need to have liminf ||pn|u n | 2 ||.L°° < IIpM 2 |U°° : to prove 
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sup {pc + pv ■ m}. 



(p\v\) 2 

this, notice that F(p,v) = p\v\ 2 /2 = is a convex functional in (pv,p) since it is 

2p 

expressed as: 

(pH) 2 

2p c+|m| 2 /2<0 

Then since p n v n = c n r] n * {pv),p n = c n {- + r] n * p) we get that 

n 

\v\ 2 

F(p n ,p n V n ) < C n 7] n * F(p,pv) < C n \\p— \\ L oo 

and letting n — > oo: 

\\d t V<f>\\ Mm < (\\p\v\ 2 \\ LO o)*C(R,d,n) 

1 

^ l|Pll|°c(B fl )ll 1; IU 00 (Br,rfp)C(- R ; ^ 

Since we impose / •<*,*) nz . 0, and sinee (1 is convex, (note that since **„ £ <\ 

a condition of this type is necessary, see [Hj, chap. 7) by Sobolev imbeddings we get also 
a bound on ||9 t $ n ||ii*(Q). This proves the first point of Theorem 12.11 

Then we obtain points 2,4,5 by the following interpolation lemma: 

Lemma 5.2 Let 4>i and <fr 2 be two R — Lipschitz convex functions on Q convex. Then 

1- there exists C,j3 > depending on (fl,R,d,p) such that 

||*i - *2lU«(n) < C||*i - $ 2||?p (n) - 

2- If moreover <&i G C 1 '" /or some < a < 1 £/ien £7iere exists C, f3' > depending 
also on a, ||$i||cn,a ; suc/i i/iai, &/ fi<s = {x G fi, ci(x, 9fi) > <5} ; 5 going to with 
||*i - *2||i,p(n), ^en 

||V*i - V* 2 |U« ( n a ) < C"||*i - *2||2, ( n)- 



Proof: Suppose that / |$i — $2| p < e p - Choose a point inside Q (say 0) such that 

Jn 

|$i(0) — *2(0)| = M. $1 and <fr 2 are globally Lipschitz with Lipschitz constant bounded 
by R. On B m / 2 r(x) HQ we have |«& a - $2 1(a) > M/2 and thus 

/ |*i - * 2 | p > vol(fi n B M/2R (x))(M/2y. 

Next note that for convex, M small enough, for any x G fi, vol(fi R Bm/ 2 r(x)) > 
Cqyo\(B m /2r(x)). Finally we have 

(33) e p > [ \& 1 -$ 2 \ p >C{n,R,d)M p+d , 

Jn 
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and thus 



M < C'(tt,R, d) 



_J B T 



which gives the first part of the lemma, with (3 



$ 1 - $ 2 | 

P 



i p 

P P+d 



p + d 



Now suppose that |V*i(0) - V$ 2 (0)| = M. One can also set $i(0) = 0, V*i(0) = 0. 
We know that $i is C 1 '" thus $i(x) < C\x\ 1+a . It follows that going in the direction of 
V3>2 one will have 



Therefore we have 



$ 2 (x) - $i(x) > M|x| - C\x\ 1+a + * 2 (0) 

r|-C|a;| 1+( 
and is equal to 



Keeping in mind that — <fr 2 (a;)| < C* e/3 yields M\x\ —C\x\ 1+a < CV 3 . The maximum 

l/a f Ar \l/a 

of the left hand side is attained for \x I 



{l+a)C 

M < Ce 13 ' 

in Qs with 5 = 6(e) going to as e goes to and with (5' 



M 

a)C 



l+a 



a/3 
l+a' 



M. 



□ 



Remark: Suppose, as it is the case for *f/, that we only know that / p|^i — ^2\ p < £ p , 
then we have instead of (f3*3*|) . 



t?> f p\^ 1 -^ 2 \P>p{B M/2R {x))MP+ d . 



The first part of the lemma yields immediately that $ G C a (I, C°(Q)) for some a > 0. 
Moreover if 0*, 02 are the Legendre transform of 0i, 2 , then ||0* — </> 2 ||l°° < ||0i — 02||l°°, 
thus * G C a (I,C°(B R )), and this gives the point 2. 

The second point of the lemma will be used to prove point 4 and 5: Indeed, if p supported 
in for some open set Q', and there exists < A, A such that A < p < A in Q', from 
Theorem 13.21 we get that for any u' CC Q', *Sf(t,-) G C 1,ai (u') for some ol\ > 0. Since 
d t $ G L U (Q), using fl2H) we get that 

J p n \dt¥ n f* < C 

uniformly in n, and thus that 

Therefore we can use Lemma 15.21 to obtain that for any uj' CC Q,\ V\£ G C^(I,C°(w')) 
(point 4 of Theorem 12 .lj) . 

Under the additional assumption that Q' is convex, Theorem 13.21 yields that <]>(£, •) in 
C 1,a2 (f2) for some a 2 > 0. The same procedure as above yields point 5. 



Now we prove the point 3 of Theorem 12.11 
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Proof of the bound on d t g 

Recall from Theorem T 



y^p„|d t V* n | < C{d,R)\\p n \\l^ {BR) \\p n vl\\l 00{BR 



We have g(t, a) = V^f(t, X(t, a)) and thus formally 

d t g(t, a) = d t W(t, X(t, a)) + D 2 V(t, X(t, a))d t X(t, a). 

Since p n converges strongly (actually weakly would be enough) to p, we know that V^ n 
converges almost everywhere to V^. (See [H] for a proof of this fact, which relies on the 
convexity of and on the uniqueness of the polar factorization). Now consider 

g n (t,a)= [ V¥ n (t,y)r tn {y-X{t,a))dy=(r in *V* n ){t,X(t,a)) 

JR d 

with i] n a smoothing kernel as above. Then g n converges almost everywhere to g. For 
/ eC°(I xfi, R d ) } let us compute 



d t g n (t, a) ■ f(t, a) dtda = 7\ + T 2 , 



with 



7i = / / / r} n {y-X{t,a))d t V* n {t,y)-f{t,a)dydadt 
J i Jn JR d 

T 2 = - [ [ [ V* n (t,y)-f(t,a)d t X(T,a)-Vri n {y-X(t,a))dydadt 
J i Jn JM. d 



Let us evaluate T\ and T 2 . 



|Ti| < / '\\f{t, .)|U~ ( n) / p{x)r) n {y-x)\d t V* n {t,y)\dxdydt 

Jl JR d xR d 

< \\f(t,-)\\L~(n)dn p n (y)\d t W n (t,y)\dxdy 

Jl JR d xR d 

< y ||/(*, •)IU 00 (^)C(-^)^)IIPn||L°°(/xlR d )ll' u n||L °(/x]R d ) 

with <i n = l/c n and from Theorem 12.21 For T 2 we have: 

mV*„(t, y) • /(*, a) d t X{T, a) ■ Vrin(y - X(t, a)) dydad< 

m9 t X'(T, a) • (,D 2 * n * r) n ){t, X(t, a)) ■ f{t, a) dydadt 
< [\\f{t,.)d t X{t,.)\\ L ^) I p{t,x){\D 2 * n \*r) n ){x)dxdt 



|T 2 



< J ||/ (t, .)d t X(t, .)|U«- ( n) ||p(*, .)IU»(R-)C(i2, d, fi) eft 

where we have used the bound on ||£) 2 ^|| L i ( Lemma l3.8|) ; we conclude that 

\\dtg\\L°°(i,M(n)) — C(R,d,£i)\\p\\L<x>(i X B R )\\dtX-\\L<=°(ixa.)- 
This achieves the proof of Theorem 12.11 



□ 
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6 Proof of Theorem 12.2 
6.1 Holder regularity 

It has been established (([13)1 and Theorem 13. lj) that d t ^f n satisfies 

V • (M„<9 t V*„) =J2 M n,iAA^n = -dtPn = ~V • (p„V n ) 

where M n is the comatrix of -D 2 \Ev To establish the Holder regularity of d t ^ n we need 
to combine three preliminary results: 

The first one (Theorem 13. 4jl asserts the Harnack inequality for solutions of the homo- 
geneous linearized Monge- Ampere equation under a condition which is satisfied when the 
density p is between two positive constants. 

The second one (Theorem 13.51 Theorem 13.71 and Corollary 13. 6j) is a local maximum 
principle that generalizes the local maximum principle for uniformly elliptic equations, to 
degenerate elliptic equations of the form V • (MVu) = V • /. The uniform ellipticity is 
relaxed to the condition that the (positive symmetric matrix valued) functions M, M~ l 
belong to LP for p large enough, p depends only on the dimension d. 

The third one ( Theorem I3.3J) asserts that the comatrix of D 2 ^/, and its inverse, are 
indeed in L^JQ!) provided that the density p is close enough to a continuous positive 
function, the closeness being measured in L°° norm. 

The result will be a consequence of the following propositions: 

Proposition 6.1 Let p = X#<ia be supported in Cl' , A and A be two positive constants 
such that < A < p(t,x) < A for all (t,x) G I xfl'. Let p n ,v n be constructed from X as 
above. Let (<& n , be associated to (p„,f2) through (00). Then there exists /3 < 1, and 
for any p > d, there exists C such that for any to' CC Q! , for any (y, r) with B ir {y) C u/, 

osc 9t * n (r/2) < (3osc dt q, n (r) + Cr s 

forn large enough. (3 < 1 depends on (A, A) (see Theorem \3.1$ , C depend on (p, A, A, inf x& / d(x, dQ'), 
\\D 2i ^\\Lp{B r { y ))), 8 = 1 - -, and 

osc u (r) = maxii — minw. 

Remark: The requirement n large enough is just to enforce that A < p n < A. 



Proposition 6.2 Under the assumptions of Proposition^^ we have, for every uj' CC Q' , 
\\dt& n \\L°°(ui') < C(K,p, inf d(x,dTl'),\,A, \\v n \\ LO o, dp u^) 

where K = ||L> 2 * n + D 2 ^- 1 ^^'), P > d. 
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Proposition 6.3 Under the assumptions of Proposition RO, for any p < oo, there exists 
e > such that if \p — 1| < e in Q' , then for every K' C Q' , K' compact, there exists Ck 1 
such that 

limsup \\D 2 * n + D 2 ^\\ LV(KI) < C K >. 



Temporarily admitting these propositions we obtain the following: 
Proof of the first part of Theorem 12.21 

From Propositions 16.11 16.21 16.31 we obtain that for any uj' CC 0! , there exists C w ', (3 < 1 
independent of n such that, for n large enough, for any B r = B r (y) C uj', with B ir C Q! , 
we have: 

osc 9t * n (r/2) < [3osc dt ^ n (r) + C w >r 5 . 

Moreover from Proposition 16.21 d t ^* n is uniformly bounded for the sup norm inside uj'. It 
is well known that this property implies Holder continuity: using ^3], Lemma 8.23, we 
obtain that for n large enough, for any w' CC Q f , there exists a > 0, C w > that do not 
depend on n, such that for any (x, y) G uj', 

\d t V n {y)-d t V n (x)\<C u ,\x-y\ a . 

Thus we have a uniform L°°(J; C a (c(j')) bound that will pass to the limit as n — > oo. We 
thus obtain the C a estimate of Theorem 12.21 

To obtain Holder continuity for d t <&, in the case where Q' is convex, we just have to use 
the identity {THD 

and the Holder regularity of V$, under the condition < A < p < A, Q' convex (Theorem 
13.2)1 . to conclude Holder regularity for 

□ 

In the next proofs we drop the suffix n for simplicity. 

Proof of Proposition fTPI This proposition is a direct consequence of Theorem 13.71 It has 
been established that d t ^ satisfies 

V ■ (MdtVV) = M iAA* = -dtp = -V • (pv) 

where M is the comatrix of D 2 *& ', given by M = det D 2 *[L> 2 *] -1 or M = pD 2 $(V*). 
We remember that < A < p < A. From Theorem 12.11 we have the a priori bound 

/ \dtQl 1 * <C(\\p n \v n \ 2 \\ L ~,n,R,d). 
Jn 
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Using then that d t ^ = — (9 t ^(V$) we have 

J P \dM u = JjdM 



and thus 



n' 



\d t *\ u <j. 



We can therefore apply Theorem 1.171 with = 1*. 

□ 



Proof of Proposition 16'. 11 

We consider a ball B ir (y) contained in Q and write d t ^ = u + w where u satisfies 



V-(MVu) = -V-H, 
u — on dB r (y), 



and w satisfies 



V • (MVw) = 

u; = dt*$f on dB r (y). 

Note that w satisfies also M^dijW = which is the equation treated in 
We denote osc/(r) = sup / — inf / and oscf(dB r ) = sup / — inf /. 

B r B r ' dB r dB r 

The assumptions of Theorem 13.41 are satisfied: indeed, in u' CC Q r , we have, for n large 
enough, A < p n < A. From Theorem 13 A\ there exists (3 < 1 such that 

osc w (r/2) < ^osCu,(r). 

From Corollary 13.61 we have 

sup |w| < C||/w ||Loor a , 

where a = 1 — d/p, C = Co||M -1 ||.ej> = Co\\p~ 1 D 2i i!\\ L p( Br ) (note that we have < A < 
p < A). Combining the two estimates, we have 

osca t v{r/2) < osc w (r/2) + osc u (r/2) 

< (3osc w (r) + Cr a 

< (3osc w (dB r ) + Cr a 

< Posc&vidBr) + Cr a 

< Posc dt *{r) + Cr a 

where in the third line we have used the maximum principle to say that osc w (r) = 
osc w (dB r ) since w can not have interior extrema. Finally we conclude 

oscQ t ^(r /2) < /3osco t iSf(r) + Cr a . 
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This achieves the proof of Proposition 16.11 □ 

Proof of Proposition \6'.cH We show here how to use the W 2,p regularity Theorem 13.31 to 
obtain estimates. First let us notice that if V\l/ satisfies (JZJ) for p supported in Q', 
< A < p < A, and since Q is convex, we know from [Sj that ^ is strictly convex in Q' 
and solution in the viscosity sense to 

det L> 2 * = p 

in Q' . Moreover ^ is in Q (Theorem 13.2)1 . From the strict convexity, for any x € Q', 
there exists a section 

S txiX = {y : *(y) < ¥(x) + V* (z).(j/ - ar) + ij 

with non-empty interior and compactly contained in Q'. (Indeed the strict convexity 
means that diameter of the sections decreases to as the height of the section t x goes to 
0). Then for every compact set K contained in Q' there exists a finite covering of K by 
sets -^Si, Si = St x ., Xi -i an d 35*% means a contraction of Si with respect to x^. Then the 
functions Ui(y) = ^(y) — U — V^ r (x i ) ■ (y — Xi) are solutions of 

det D 2 Ui = p in Si 
Ui = on dSi. 

From John's lemma (see [Hj), we can find an affine transformation Tj, with detTj = 1 
and a real number /ij such that B\ C p^ T^ x [Si) = Si C dBi. Finally, considering 
Wi(y) = \ui{p{Ti y) we get that is solution to 

detD 2 £ti(y) = p(y) = p(piTiy) in ^ 

Mj = on dS; 

Si(xi) C5(C dB^Xi). 

We can invoke Theorem 13.31 for For any < p < oo, if \p — 1| < e(p) (this property is 
invariant under the renormalizations performed above), we have 

\\D 2 Ui + D 2 uT l \\ LP[Bl) <C 

(measO^r 1 ^ \\D 2 Ui + DVlLn^s,) < C\\Ti\\ 2 . 

By our covering process, we have K c{J { TiPiBi (xj). It follows that for every compact set 
K C Q', there exists and constant Ck such that HD 2 ^^^) < Ck and || Z? 2 ^ - 1 1| ^(.fs") < 
Ck- The constant depends on the supremum of the norm of the transformations T; 
and can be taken (by compactness) uniformly bounded given Q, Q', K, A, A. 

Now we show that this covering process behaves uniformly well when we consider the 
regularization p n of p and let n go to oo. Indeed the corresponding \I> n will converge 
uniformly to ^ and since the limit \I/ n is C 1 the sequence V\I/ n converges also uniformly 
in every compact set of Q'. Therefore the set S 1 " = {y, ^f n (y) < ^n(^i) + V^nO^i) ■ (y — 
Xi) + ti} converge uniformly to Si. This means that for n large enough, the set K will be 
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covered by [j { ^5*™. Consider /x™, T™ the corresponding normalization, then we also have 
T™,/!™ converging to T i: p^, and K will be covered by ^T^p^Bi (xi). 

Moreover since we consider a compact set K contained in Q' and since \p — 1| < e in 
O', it follows from the construction of p n that, for n large enough, \p n — 1| < e in For 
n large enough, the functions (obtained by the renormalization procedure) will thus 
all satisfy the assumptions of Theorem 13.31 

Therefore, for every K CC Q', there exists Ck independent of n such that, for n large 
enough, 

\\D 2 ^ n + D 2 ^- 1 \\ LP{K) < C K . 
This achieves the proof of Proposition 16.31 

□ 

Proof of the gradient bounds 

This is point 2 of Theorem 12.21 The gradient bounds follow directly from Proposition 15.11 
combined with Proposition 16 .31 In estimates ()27|l2T)|) take r = oo. Note that from Lemma 
13.81 we have the bound HD 2 ^^^) < C(R,d,fl). This ends the proof of Theorem 12.21 

□ 

7 The periodic case: proof of Theorem 12.3 

This result is only an adaptation of the two previous Theorems. All the regularity results 
used adapt to the periodic case as follows: 

Theorem 7.1 Let p be a Lebesgue integrable probability measure on M d /Z d . There exists 
a unique convex on M. d , with — \x\ 2 /2 periodic, that satisfies 

det £> 2 * = p 

in the following sense: 

V/ G C°(R d /Z d ), [ p/(V*) = / /. 

It has the following regularity properties: 

1. If for some pair (A, A) G M^_, we have A < p < A, then for some a > depending 
on A/A, * - |x| 2 /2 is m C^T*). 

2. For every p < oo, there exists e(p) such that if \p — 1| < e(p) ; then ^ — |x| 2 /2 G 
W 2 ' p (T d ). 

3. If p is positive and in C°°(T d ), then * - |x| 2 /2 G C°°(T d ). 
We then modify the approximation procedure as follows: we take 

* = c,<*., + ±) 
n 

det£> 2 * n = p n 
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with the constant c n such that j Td p n = 1. Then we use the same techniques as in the 
Theorems OO 

We only mention the two new results that arise in this case: 

In point 2, we obtain that g G C a (I, L°°(T d )). Indeed, g = V*(t,X(i)). We al- 
ready know that, under the present assumptions, G C a (I x T d ), moreover X G 
W l '°°(I : L°°(T d )) and the result follows. 

In point 4, under the assumption that \\p— 1 ||i°°(/xT d ) < e for e small enough depending 
on q, , we are able to obtain a bound in L q (T d ), q < 2 for d t g. Indeed, writing 

g n {t,a) =V* n {t,X{t,a)) 

as in the proof of Theorem 12.11 and differentiating with respect to time, we obtain 

d t g n (t, a) = d t VV n (t, X(t, a)) + D 2 * n (t, X(t, a))d t X(t, a). 

with ^ n obtained from p n , and thus in C°°(I x T d ). If p is close enough to 1 so that 
D 2 ty n is bounded in L p (T d ) (cf. Theorem 17. II above) . the first term is bounded in L q (T d ), 
with q = (as in Proposition 14. ljl . The second term is bounded in L p (T d ). Then we 
let g n converge to g. 

Note that this bound can not be obtained in the non periodic case since we have only 
interior regularity available for 

8 Counter-examples 

Here we show through some examples that the bounds obtained in Theorem 12. II are sharp 
under our present assumptions. 

Example 1: d t V® (£ L} oc and d t $> £ C°. 
Consider in Q = B(0, 1) in IR 2 , and X(t, •) : B(0, 1) — > IR 2 defined with complex notations 
X = x + iy by 

on y > 0, 

X(t, (x,y)) = e it (x + ty)+it, 

on y < 0, 

X(t,(x,y)) = e it (x + iy)+t 2 . 

We check that X#rfa has a density bounded by 1, that d t X G x R+). If X = V$ 051 

is the polar factorization of X then up to a constant, $ is defined for t > 0, (x, y) G Q by: 

*(t, (x, y)) = sup{i(x 2 + y 2 ) + t 2 x, i(x 2 + y 2 ) + ty}. 
On {y > tx} we have 



$(t,(x,y)) = -(x 2 + y 2 )+ty, 
V*(t, = (x,y) + (0,t), 



8 COUNTER-EXAMPLES 



28 



and on {y < tx} 

&(t,(x,y)) = l -{x 2 + y 2 )+t 2 x, 
V*(t, (x,y)) = (x,y) + (t 2 ,0). 

Thus 

d t ®(t, (x,y)) = yx{ y >tx} + 2txx{ y< tx} i C°, 

ftV*(t, (x,y)) = (0, l) X {y>t x} + (2t,Q) X{ yM + (t 2 , -t)H d - l {y = tx} $ L} oc 

Example 2: Here we adapt a counterexample of Wang to build an example of a 
solution where dt^ ^ C°. 
In R d , let x = (a;j)i<j<d and 

X(0,ac) = V$ (x) 

$>o( x ) convex Lipschitz on fl, <fr = +oo outside, such that p = V&o(x)#dx has a density 
in L°°(R 2 ). Let 

X(t,x) = V®o(x) + tv 

for some fixed v G IR d . X is Lipschitz with respect to time. Then 

<&(t, x) = &(x) + ta; • v, 
S/$(t,x) = V$ (» + 

If is the Legendre transform of $o, the Legendre transform of •) is given by 

= * (z - tv), 
V*(t,x) = V* (^-to), 

thus 

d t ^{t, x) = v ■ V^ {x - tv), 
d t W(t, x) = D 2 & (x - tv) ■ v. 

Wang has shown in j2H] some counterexamples to the regularity of solutions of Monge- 
Ampere equations: namely, for d > 3 he has exhibited a solution u of 

det D 2 u = f 

with / only bounded by above, such that u C 1 . By taking SE^ = u one has an example 
of time dependent map such that 

d t ^{t, x)=v W (x - tv) (£ C°. 
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9 Application: the semi-geostrophic equations 

The semi-geostrophic system is derived as an approximation to the primitive equations in 
meteorology, and is believed to model frontogenesis (see [12] )■ The formulation of the 3-d 
incompressible version is the following: we look for a time dependent probability measure 
p that satisfies the following SG system: 



Here v 1 - means (— v 2 ,vi,0). Equation (j3Tlj) is understood in the sense of (JJJ), where an 
open set ft of total mass 1 has been given before. 

The system has also a periodic version in which ft = T 3 itself and equation (|37)|) is solved 
with the condition that ^ — \x\ 2 /2 is Z 3 periodic. 

The set ft is here called the physical space, whereas the space in which p lives is the dual 
space. Existence of global weak solutions for the SG system with initial data in L 1 has 
been proved in [2|, JI] and JHj- Note that uniqueness of weak solutions is still an open 
question. 

9.1 The Lagrangian formulation of the (SG) system 
Here we look for a mapping X : R + x ft -> 1R 3 that satisfies 



If we define p(t) = ~X(t)#da, the last equation means that for all t, Vf'(t) solves det D 2 ^f(t) = 
p(t) in the sense of (JTj). Having X solution of (|3*71 l3*Sjl implies that p(t) = X(t)#da is 
solution of (|3*3l 1331 l3Tlj) . X defines the characteristics in the dual space whereas g defines 
the characteristics in the physical space. 

We expose briefly the arguments that allow to define the characteristics of the SG 
system: 

1- First we check that X(i) will satisfy for any time t the condition (JIJ: indeed, the flow 
being incompressible, all the LP norms of p are conserved. Therefore, given the potential 
*&{t), if X satisfies the condition (0), or equivalently if po G L 1 , then we know a priori 
that X(t) satisfies the condition JJJ for all time. 

2- The velocity field is a priori bounded in BV because of the convexity of St (see Lemma 
13.8)1 . Moreover it is incompressible. Therefore thanks to the result of pQ, the characteris- 
tics of the corresponding ODE are uniquely defined for almost every initial data, which 
means that the curve t \— > X(t, a) is uniquely defined for almost every a e Q. 

For Q bounded, it is easily checked (see [2j) that if X G L°°(Q), then (X, d t X) e 
L°°([0,T] x ft) for all T > 0. The velocity field being incompressible, if p e L°°(E 3 ), 
then p G L°°(]R + x IR 3 ). Note that the Lagrangian system can also be defined in a periodic 
space, where X is periodic in space for all time, and we require ^ — \x\ 2 /2 to be periodic. 
The bound of X, <9 t X in L°°(IR + x T 3 ) is then independent of the initial data. Moreover, 
in this setting, if p is such that 



(34) 
(35) 
(36) 



d t p + V • (pv) = 
v(t,x) = (V*(t,x) -x) 
det D 2 ^f(t,x) = p(t,x). 



± 



(37) 
(38) 



d t X(t, a) = (V*(t, X(t, a)) - X(t, a)) x 
V*(t) o X(t) = g(t) G G(ft), * convex. 



(39) 



< A < po < A 
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for two constants A, A, this property remains satisfied for all time, once again due to the 
incompressibility of the velocity field. 
Thus we conclude the following: 

Lemma 9.1 Let X G L°°(fi;R 3 ), p = X #da G L°°(R 3 ). Then p, X the corresponding 
solution of the SG system satisfies for all T > 0, 

X,<9 t X G L°°([0,T] x O) 
p G L°°(R + x R 3 ). 

In the periodic case this remains true, and if moreover po satisfies then for all time 
t, p{t) satisfies < f3^) . 

Under the assumptions of the above lemma, it is clear that X satisfies the assumptions of 
Theorem 12.11 In the periodic case, if satisfied at time 0, all the assumptions of Theorem 
12 .31 are satisfied for all time. We can now state the following theorem of partial regularity. 
We restrict ourselves to the periodic case. 

Remark: We also conjecture that the assumptions of Theorem 12 . 21 can be satisfied for some 
finite time, but the control the evolution of the support of p poses some some difficulties. 

Theorem 9.2 Let X, p, g, <fr be as above, with p = X.#da be a space-periodic solution 
°f EH \HB\) , an d X the corresponding space-periodic solution of 135)) . Suppose that 
p G L°°(T 3 ), then 

^gGL°°(M + ,.M(T 3 )), 
d t V$ G L°°(R + ,.M(T 3 )). 

// moreover there exists < A, A such that \ < Po < A, then there exists a > depending 
on (A, A) such that 

g G C a (R + ,L°°(T 3 )). 
For all p < 2, there exists e(p), such that if |p — 1| < e G°); then 

^gGL-([0,T],Z?(T 3 )). 
There exists e , such that if \p — 1| < e < e 0; then 

d t <f>,d t V G L°°(R + ,C Q (T 3 )) 

where a > depends on e. 

Remark: The equations of motion in physical space We derive here formally 
the equation giving the evolution of g: writing (|T0|) with v as above, we have 

(x - V*) 1 " = v(V$) = d t V<$> + D 2 <f>w : 
V- w = 0, 
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where 9 t g(g _1 ) = w. This equation formally determines the evolution of the system, since 
the knowledge of 3>(i) determines a unique pair 9 t V$, w satisfying the above decomposi- 
tion (see Proposition II. 1|) . One can see a parallel with the Euler incompressible equation 
where the evolution is given by solving the following decomposition problem: 

— v ■ Vf = d t v + Vp, 
V • v = 0. 

Thus the semi-geostrophic equations are associated to the decomposition of vector fields 
of Proposition 11.11 in a similar way as the Euler incompressible equations are associated 
to the Hodge "div-curl" decomposition. 
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